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We present a mass spectrometric technique for identifying the masses and relative 
abundances of Coulomb-crystallized ions held in a linear Paul trap. A digital ra- 
diofrequency waveform is employed to generate the trapping potential, as this can be 
cleanly switched off, and static dipolar fields subsequently applied to the trap elec¬ 
trodes for ion ejection. Excellent detection efficiency is demonstrated for Ca + and 
CaF + ions from bi-component Ca + /CaF + Coulomb crystals prepared by reaction of 
Ca + with CH 3 F. A quantitative linear relationship is observed between ion num¬ 
ber and the corresponding integrated TOF peak, independent of the ionic species. 
The technique is applicable to a diverse range of multi-component Coulomb crystals 
- demonstrated here for Ca^/Nldg/NH^ and Ca + /CaOH + /CaOD + crystals - and 
will facilitate the measurement of ion-molecule reaction rates and branching ratios in 
complicated reaction systems. 
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I. INTRODUCTION 


Recent improvements in cold molecular beam methods and their combination with ion 
traps is facilitating exciting new developments in the observation of cold, controlled ion- 
molecule reactions. 12 It is not sufficient, however, to simply conduct experiments in the cold 
regime; the resulting products must be quantitatively detected to extract information about 
the reaction process, such as the rate, product branching ratio, and energy dependence. 

Trapped, laser-cooled ions can undergo a phase transition to adopt a highly ordered 
‘Coulomb crystal’ structure. Co-trapped ions of a different species may be sympatheti¬ 
cally cooled through elastic collisions with the laser-cooled ions, forming multi-component 
Coulomb crystals. In this way, cold molecular ion targets of species that cannot be laser 
cooled can be prepared. Imaging the fluorescence continuously emitted by the laser-cooled 
ions enables real-time observation of the crystal framework, as depicted in figure [jjb). Due 
to the dependence of the trapping held on the charge-to-mass ratio of the ions, different 
ionic species are located in separate regions of the Coulomb crystal. In chemical reac¬ 
tions, this separation can be exploited and reaction rates measured through monitoring the 
changing shape of the crystal. For instance, the reaction of laser-cooled ions with neutral 
reactants - and the disappearance of fluorescing ions - is accompanied by changes in the 
crystal framework as product ions heavier (lighter) than the laser-cooled species migrate 
to the outer (inner) region of the crystal . 13 Hence, even though the non-fluorescing sym¬ 
pathetically cooled ions can not be observed directly, comparison of molecular dynamics 
(MD) simulations with observed fluorescence images enables quantitative determination of 
ion numbers (see Fig. [Tj). 

In complicated chemical systems with multiple product channels, however, MD simula¬ 
tions alone are insufficient for establishing accurate reaction rates. For example, in the re¬ 
action between C 2 H^ (sympathetically cooled into a Ca + Coulomb crystal) and NH 3 , there 
are numerous exothermic pathways: proton transfer, charge transfer, two H-elimination 
pathways and two ^-elimination pathways.® In this case, both the reactant ions and nu¬ 
merous ‘dark’ (non-fluorescing) product ions are lighter than the laser-cooled species, and 
thus would be indistinguishable in the dark inner core. A change in the ratio of the num¬ 
bers of light ions of different mass would have little effect on the crystal framework. In 
some instances, resonant-excitation mass spectrometry (MS) methods can be utilized. The 
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FIG. 1. The number of ions in a given Ca + Coulomb crystal can be established to within ±5 ions 
by fitting the experimental image to MD simulations. A simulated crystal image with 50 Ca + ions 
is shown in panel (a), alongside an experimental Ca + crystal image in panel (b). The simulated 
crystal is overlaid on the experimental image in (c), illustrating the excellent agreement between 
the positions of the ions in the simulated and experimental crystals. This approach is also adopted 
for multi-component crystals; in panel (d), a bi-component Ca + /CaF + experimental crystal (in 
which only the Ca + are fluorescing) is overlaid with a simulated Ca + /CaF + crystal comprising 
30 Ca + and 35 CaF + ions. Only the Ca + ions are made visible in the simulation, for ease of 
comparison between the images. 


scanning of an additional radiofrequency (RF) field resonantly excites the secular motional 
frequencies of trapped ions, allowing the ion masses to be determined. 5 ' In the limiting 
case of a single molecular ion sympathetically cooled by a single laser-cooled ion, the mass 
of the molecular ion can be precisely measured, with for example 26 MgH + and 26 MgD + ex¬ 
perimentally distinguishable. 8 While several variants of this non-destructive technique have 
been developed]- 9 10 resonant-excitation MS remains limited in application to small crys¬ 
tals with only a few different ion types; when motional frequency spectra are recorded for 
large multi-component crystals, the features are poorly resolved, difficult to unambiguously 
assign, and the relative abundances cannot be ascertained quantitatively)- 3 An alternative 
approach, termed ion crystal weighing, employs a pulsed electric field to excite the center- 
of-mass mode of the ions. 11 The frequency of the mode, and hence the mass of the crystal, 
can be deduced from the Fourier-transformed autocorrelation of the ion fluorescence signal. 
While ion crystal weighing is applicable to larger crystals, interpretation of the results is 
challenging when more than 2 ionic species are present or when the identities of the dark 
ions are unknown. 

Adopting the principles of Wiley-McLaren time-of-flight MS (TOF-MS)) 12 Hudson and 
co-workers recently demonstrated the ejection of ions from a linear Paul trap into a TOF tube 
and onto an external detector. 13 Ion ejection was achieved by switching off the RF voltages 
and converting the trap electrodes into a repcllcr-extractor pair by applying static voltages. 
However, when using resonantly-driven cosine RF voltages, the trap electrodes form part of 
a resonant inductor-capacitor circuit. As the voltage is switched from cosine RF to direct 
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current (DC), the amplitude of the trapping voltages decays with a time constant 1/k on 
the order of microseconds, and ringing of the form V r fe~ kt cos(iut) is superimposed on the 
ejection voltages. This compromises detection efficiency; ions experience vertical deflection 
away from the TOF axis and typically miss the detector. In an attempt to overcome this 
problem, Hudson and co-workers applied very high repcller and extractor voltages (1.4 kV 
and 1.2 kV, respectively). The observed mass-spectrum peak intensities, however, were not 
representative of the relative quantities of each species.^ 


Quadrupole ion traps have been successfully combined with mass analysis apparatus in a 
variety of configurations - with ion ejection both along and perpendicular to the ion trap axis 
- and many such instruments are now commercially available. 14 15 Such systems, however, are 
not designed for the study of ion-molecule reactions in Coulomb crystals. There is insufficient 
optical access to the trap for effective laser cooling and crystal imaging. Furthermore, the 
typical operating pressures and repeated ejection cycles are not appropriate for the long 
storage times required in our experiments - we require ultra-high vacuum conditions to 
achieve ion storage times on the order of tens of minutes. 


In this paper, we present a new digital ion trap (DIT) Wiley-McLaren MS technique 
for identifying the masses and relative abundances of all ions within a Coulomb crystal. 
Digital trapping waveforms, also termed square-wave or pulsed waveforms, are applied to 
a linear Paul trap in place of the conventional cosine waveform. This facilitates rapid, 
clean switching between the RF trapping and static DC ejection voltages, with no residual 
ringing. The entire crystal is ejected radially from the trap, where the ions pass through 
a TOF tube onto a microchannel plate (MCP) detector for mass-sensitive detection. The 
quantitative performance of the technique is characterized by the analysis of hundreds of Ca + 
and Ca + /CaF + Coulomb crystals, in addition to multi-component Ca + /CaOH + /CaOD + 
and NHg/NH 4 /Ca + crystals. The efficiency of ion ejection is studied, aided by extensive 
modeling. The methodology introduced in this work is applicable to any reaction of laser- 
cooled or sympathetically-cooled ions, provided two conditions are satisfied: the co-trapped 
ions have an appropriate mass-to-charge ratio for efficient sympathetic cooling; and any 
product ions are formed with insufficient kinetic energy to escape the trap. 
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II. EXPERIMENTAL METHODS 


The experiments presented in this paper utilize a linear Paul trap operated in ultrahigh 
vacuum (UHV) conditions, illustrated in Fig. [2j The cylindrical trap electrodes, radius 
4 mm, have a diagonal electrode-surface separation of 2ro = 7.0 mm and are each composed 
of three segments (endcap separation 2zo = 5.5 mm), to which a combination of radiofre¬ 
quency (RF) and static voltages are applied. The conventional time-dependent trapping 
potential in a linear Paul trap is of the form cp r{ (x,y,t) = ~ 2 y j cos(u; r ft), with l/ rf the 

peak-to-peak voltage amplitude, ca r f the RF drive frequency. However, as first demonstrated 
in the 1970s]a cosine wavefunction is not essential for the generation of a trapping po¬ 
tential. A DIT, as we utilize in this work, employs a rectangular waveform of period T and 
fractional pulse width r, yielding a time-dependent trapping potential 18 19 
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where we can represent the pulsed waveform as the sum of its cosine Fourier components, 

OO 

P T (f) = a n cos(?ro; r f*). 

n= 1 

Axial confinement of ions in the trap is obtained through the application of static voltages 
(P dc ) to the endpieces of the segmented electrodes, <f) en d(x,y, z) = (z 2 — - ~t, y j, with 

t] = 0.244 a geometrical factor for this trap. The total electric potential is thus </> = 
0 r f(x, y , *) + 0end(A> V, z). The equations of motion for a single trapped ion can be expressed 
as a Mathieu equation, d 2 u/d£ 2 + u[a u + 2q u f(C,)\ = 0, where /(£) = cos(2£) for a cosine 
RF field and /(£) = P r (0 for a digital RF field, u E [x,y] and £ = ^ca r f*. The stability 
of an ion trajectory in the trap is dependent on the dimensionless Mathieu parameters a u 
and q u . The ( q u , a u ) plane is divided into stable and unstable regions; for stable values of 
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a u and q Ul the solutions to the Mathieu equations are bound and the ion will be confined 
in the trap. In a DIT, the Mathieu equations are dependent on the fractional pulse width 
(r), which can be varied to yield a stability diagram comparable to that of a cosine trap. 
(See reference p3] for a derivation of the stability criterion in a DIT.) Bandelow et al. have 
recently experimentally confirmed the applicability of the a u and q u parameters in mapping 
the stability regions of a hyperbolic Paul DIT. 20 



FIG. 2. (color online). Schematic illustration of the linear Paul trap utilized in this work, with four 
electrodes in a quadrupole arrangement. Each electrode is composed of three segments, with RF 
voltages ( V T {) applied to all segments and additional DC voltages (U dc ) applied to the endpieces. 


We can define dimensionless Mathieu parameters which are dependent on the mass (m) 
and charge ( Q ) of the species, 


&X Q"y 2 
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(4) 


q x = -q y = - 2 ^ V, 'i ; q z = 0, (5) 

which enable one to construct a DIT stability diagram in terms of U dc and Vfi (see Fig. [3]). 
From the secular potential, defined as the sum of the pseudo-potential (the harmonic time- 
averaged potential arising from the oscillating RF fields) and the static endpiece potential, 
one can derive the digital trap depth, 
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OO 

with c=;E a n/ n2 a sum over the Fourier-component amplitudes a n (in a cosine trap, 

71=1 

c = 0.5). The digital RF voltage is applied at a frequency of oj t f = 27r x 1.33 MHz, with 
typical trap parameters F)f = 50 V, U dc = 1.0 V and r = 0.25 that give Matthieu parameters 
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falling within the regions of stability for all ionic species considered in this work. 

The experimental digital waveform (Fig. [4]) exhibits a non-zero pulse decay time, and thus 
deviates from the ideal digital waveform. MD simulations - in which the measured digital 
waveform can be directly utilized - indicate that this deviation from the ideal waveform 
has a very minor impact on the trap depth and does not adversely affect crystal formation, 
crystal stability or the ejection properties of the system. This is in agreement with Sudakov 
et a/., who found non-ideal waveforms to have an insignificant effect on the calculated regions 
of stability in a DIT. 21 Hence assuming the real digital waveform is well represented by an 
ideal digital waveform, a trap depth of ~1.21 eV is calculated - shallower than the cosine 
trap depth of ~1.36 eV as previously utilized in our research group^, but entirely sufficient 
for most applications. 



0 100 200 300 400 0 50 100 150 

V rr (V) v rf (V) 


FIG. 3. (color online), (a) The regions of stability for several ionic species in the DIT operating 
at t = 0.25 and cj r f = 2ir x 1.33 MHz. For clarity, only the stability region of Ca“*" is shaded black; 
for all other species, the stability region is enclosed between the relevant colored (shaded) curves. 
Regions of overlapping stability indicate the operating parameters (Vj-fj Hdc) where species can be 
co-trapped. (b) A closer view of the experimentally-accessible region of the stability diagram. 



FIG. 4. Experimental digital RF trapping waveform. 

To form Coulomb crystals, calcium atoms in an effusive skimmed beam are non-resonantly 
ionized (at 355 nm) in the trap center. The resulting Ca + ions are laser cooled on the 
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4s 2 S' 1/2 —> 4 p 2 P \/2 transition at 397 nm, with a 866 nm repump beam addressing population 
lost to the 3d 2 D 3 / 2 state. The cooling lasers are reflected back along the trap axis to achieve 
bi-directional cooling. An additional 397 nm beam at 45° to the trap axis provides radial 
cooling. Fluorescence continuously emitted by laser-cooled Ca + ions is imaged by a lens into 
an intensified charge-coupled device camera system, yielding a two-dimensional image of the 
Coulomb crystal. High-precision leak valves facilitate the controlled introduction of neutral 
reactants into the ion trap chamber, typically held at pressures ~ 1 x ICC 9 mbar. In this 
way, we generate molecular ions through chemical reactions between the neutral reactant 
gas and trapped Ca + ions, or with trapped sympathetically cooled ions. 

After the quadrupole trapping voltages are switched off, ions are radially ejected into the 
TOF tube by dipolar repellcr and extractor fields of +220 V and +97 V respectively, applied 
with a rise time of 70 ns, and are accelerated towards a grounded mesh of 88% transparency 
mounted across the entrance of the TOF tube (see Fig. |5](b)) . The ions are detected by a 
microchanncl plate (MCP) detector at the end of the flight tube. The repellcr and extractor 
voltages employed were found empirically to give the best TOF mass resolution for Ca + , 
and this two-stage acceleration arrangement is equivalent to a Wiley-McLaren time-of-flight 
scheme 12 . In the experiments reported here there was (nominally) zero delay in the switching 
from the trapping voltages to the ejection voltages; although the delay can be varied easily 
in the digital trap, minimizing the delay was found to optimize the mass resolution. 


III. SIMULATIONS 

Coulomb crystals are simulated using a custom-written MD crystal simulation program, 
which takes the experimental digital waveform as the input RF trapping voltage. Simulated 
crystal images can thus be directly compared to those observed experimentally, as demon¬ 
strated in Fig. [lj A MATLAB program is used to model the ion trajectories upon ejection 
by calculating the electric fields in the ion trap and TOF tube, using SIMION to solve the 
Laplace equation numerically and including ion-ion repulsion. Taking the identity, position 
and velocity of each ion within a Coulomb crystal (as output by the MD code) and the 
experimental ejection voltages as input parameters, ion trajectories are calculated using a 
Velocity-Verlet algorithm. For Coulomb crystals comprised of < 200 ions, simulations indi¬ 
cate that the detection efficiency is limited only by the transmittance of the grounded mesh; 



32 mm 268 mm 


FIG. 5. (color online), (a) The integrated Ca + peak intensity in the TOF spectrum is plotted for 
each Coulomb crystal against the number of Ca + ions, determined from fitting the experimental 
images to MD simulations, (b) Schematic illustration of the TOF-MS apparatus, depicting the ion 
trap electrodes (viewed along the trap axis, with ion ejection perpendicular to this axis), TOF tube 
and MCPs. The electric potential near the ion trap is positive (shaded red online), and near the 
front detector plate is negative (shaded green online) with the trajectories of 200 Ca + ions ejected 
from the trap shown (in blue online.) 


detection efficiency is independent of the particular species, for all species considered in this 
work. 


IV. RESULTS AND DISCUSSION 

To verify the predicted detection efficiency experimentally, the number of Ca + ions in a 
given crystal is determined for several hundred Coulomb crystals of various sizes by compar¬ 
ing the experimental image immediately prior to ejection to a set of simulated crystal images 
(see Fig. [I]). Following crystal ejection, the area under the relevant mass peak is calculated 
by integration of the (background-subtracted) TOF trace. Fig. |5](a) illustrates the strong 
linear correlation between the number of ions in a Ca + Coulomb crystal and the integrated 
TOF peak recorded upon crystal ejection, for up to 200-ion crystals. The measurement error 
is governed by the uncertainty in the comparison between an experimental crystal and its 
MD simulation (±5 ions). The resolution of the mass spectrum, defined as m/ Am where 
Am is the FWHM of the mass spectrum peak, is 90 for a typical 100-ion Ca + crystal and 
> 70 for all Ca + crystals containing up to 200 ions: The slightly lower resolution for larger 
crystals is due to the greater spread of initial positions and velocities in those cases. Crystals 
with more than 200 ions were not considered experimentally due to imaging limitations; the 
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objective lens and camera restrict the field of view. 

To establish the applicability of this DIT MS technique to multi-component Coulomb 
crystals, mixed crystals of Ca + and CaF + are prepared by admitting a low pressure of CH 3 F 
gas to the ion trap chamber, leading to the reaction Ca + + CH 3 F —y CaF + + CH 3 . The 
resulting CaF + product ions sympathetically cool into the Ca + Coulomb crystal, forming a 
dark outer shell - indicated by the flattening of the fluorescing Ca + core in the crystal images 
(see Figs. 0 d) and[b](a)). Repeating simulations for variable numbers of laser-cooled ions 
and dark ions enables a determination of the number of each species in the crystal, assuming 
there is only one “dark” component. The mass-to-charge ratio of CaF + :Ca + is 1.475, well 
within the desired range for optimal sympathetic cooling . 1 2 To our knowledge, this represents 
the first demonstration that multi-component Coulomb crystals can be formed in a DIT. 
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FIG. 6. (a) A bi-component Ca + /CaF + Coulomb crystal and the TOF trace recorded follow¬ 
ing ejection of the crystal. Peaks arising from species with masses differing by 1 amu can be 
resolved, as demonstrated by the ejection of (b) Ca + /CaOH + /CaOD + and (c) NHg“/Nfp|"/Ca + 
multi-component crystals. The two unlabelled small peaks in (a) are due to unidentified minor 
impurities. 


As Fig. [7] illustrates, a linear relationship is observed between the integrated TOF 
mass peaks and ion numbers - for both the fluorescing Ca + and the non-fluorescing CaF + 
ions. The gradients of the lines-of-best-fit, (4.0 ± 0.2) xlO -10 Vs ion -1 for Ca + and (3.9 ± 
0.4)xlO -10 Vs ion -1 for CaF + (for up to 60 CaF + ions), established using a linear least- 
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squares fitting procedure with errors considered in both coordinates, are in quantitative 
agreement. These results are also in excellent agreement with the previous Ca + -only crystal 
results shown in figure [5][a)) , demonstrating the utility of the technique: equal detection 
efficiency is observed regardless of ion mass. There are too few crystals in the data set with 
> 60 CaF + ions to extend the trend further, and the precision with which CaF + ion num¬ 
bers can be estimated from simulations decreases at higher CaF + :Ca + ratios. Trajectory 
simulations indicate that the observed linear trend will hold for larger numbers of ions (up 
to 400-ion crystals), and for a diverse range of ionic species, provided the mass-to-charge 
ratios of co-trapped species are similar and there are a sufficient number of laser-cooled ions 
for efficient sympathetic cooling. 




Number of Ca + ions 


Number of CaF + ions 


FIG. 7. Bi-component Ca + /CaF + crystal ejection analysis. The integrated TOF peak for (a) 
Ca + and (b) CaF + ions is plotted against the number of ions of that species in each crystal, 
ascertained from fitting the crystal images to simulations, for bi-component crystals of various 
sizes and compositions. 


The mass resolution of the technique is verified by the ejection of crystals containing 
species differing in mass by 1 amu. Tri-component Ca + /CaOH + /CaOD + crystals are formed 
in a two-step process. H 2 O vapor is admitted to the chamber through a leak valve, where 
reaction with Ca + (in the 3d 2 D excited states populated by the laser-cooling process)^ 
yields CaOH + molecular ions. D 2 O is then admitted through a second leak valve, forming 
CaOD + ions. As figure | 6 |b) illustrates, peaks arising from CaOH + and CaOD + are resolved 
in the spectrum, indicating a resolution better than 1 mass unit at 58 amu. 

The applicability of this DIT MS technique for species lighter than Ca + is confirmed with 
tri-component NH^/NHj"/Ca + Coulomb crystals. NH 3 is resonantly ionized at 312.5 nm, 
with some of the sympathetically cooled Nhhj“ ions subsequently reacting with neutral NH 3 
to form NHjJ", yielding tri-component crystals. Following the same imaging and ejection 
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process detailed above, the TOF spectra again exhibit well-resolved peaks, as shown in figure 
[6^c). In this system, both the dark species are lighter than Ca + . Accordingly, conventional 
comparisons of the experimental crystals with simulated images do not reveal the relative 
numbers of NH^ and NH^ ions, as both species are located in the dark inner core. The 
same applies to multi-component Ca + /CaOH + /CaOD + crystals, where the CaOH + and 
CaOD + species are both heavier than Ca + and therefore are both located in the dark outer 
shell. DIT MS enables us to acertain the relative numbers of dark ions which are otherwise 
indistinguishable in the experimental and simulated images. 


V. CONCLUSIONS 

In summary, we have established DIT MS to be a robust technique for the quantitative 
characterization of multi-component Coulomb crystals produced, in this work, by chemical 
reactions of the trapped ions. The ejection process is efficient, with a linear correlation 
between the number of ions ejected and the corresponding integrated TOF peak, and a 
detection efficiency that is independent of the crystal size or the ionic species involved. The 
DIT MS technique is applicable to a diverse range of multi-component Coulomb crystals 
and hence to a diverse range of chemical reactions - not only reactions of laser cooled ions 
but also of sympathetically cooled ions. While the ejection process is destructive, repeating 
the process with reloaded crystals at a range of reaction times allows the determination of 
accurate reaction rates and branching ratios in complicated ion-molecule reactions. With the 
ability to unambiguously detect both the masses and relative numbers of trapped ions, DIT 
MS has the potential to be a powerful detection technique in the study of sympathetically- 
cooled molecular ion reactions - utilized as a stand-alone method or in combination with 
existing techniques such as MD simulations. 
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